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ABSTRACT 


In  this  paper  an  attempt  has  been  made  to  clarify  certain 

[2] 
questions  about  the  convergence  of  the  Bremmer  series  solution    of 

the  linear  second  order  differential  equation 


u"  +  k^(x)u  =  0 


These  questions  have  been  brought  forth  recently  in  mathematical 

literature,  in  particular  in  an  article  by  F.  V.  Atkinson    and  in 

[51 
an  earlier  one  by  Bellman  and  Kalaba.     It  is  shown  here  that  series 

of  the  Bremmer  type,  other  than  that  discussed  by  Atkinson,  exist  and 

have  better  convergence  properties.   Moreover,  such  series  have  been 

fftl 
given  by  J.  B.  Keller  and  H.  B.  Keller    for  the  more  general  case  of 

a  system  of  first  order  differential  equations - 
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I.  Introduction 

F.  V.  Atkinson^ L J  in  an  elegant  analysis,  has  investigated  the 
validity  of  a  multiply  reflected  wave  series  solution  of  the  differential 
equation 

u"  +  k^(x)u  =  0  .  '  (1.1) 

The  solution  has  the  form 

u(x)  =  u^(x)  +  u^(x)  +  u^Cx)  +  ...   ,  (1.2) 

each  term  of  which  corresponds  to  a  wave  traveling  to  the  left  or  to  the  right. 
The  solution  is  assumed  to  be  one  such  that  the  sum  of  all  the  left  moving 
waves  vanishes  in  the  limit  x  ->  +oo.   The  terms  are  classified  according  to 
the  number  of  reflections  the  corresponding  wave  has  undergone. 

•  The  solution  (1.2)  was  given  independently  by  Bremmer^L  -I  Landauer^L  J 
and  Schelkunoff .  L  J   it  was  discussed  by  Bellman  and  Kalabai-  J  who  also 
considered  its  validity.   Atkinson!-  J  showed  that  if  k(x)  is  continuously 
differentiable  and  if  k(xj  satisfies  the  conditions 

00 

k(x)>a">  0,      /  |k'(x)  |dx  <  00  ,  (l-5-^) 

o 

where  a  is  a  constant,  as  well  as 

00 

J  |k'(x)/k(x)|dx  <  «  ,  (1.5) 

o 

the  series  (1.2)  converges.  He  also  showed  that  the  number  jt  in  (1.5)  cannot 


The  second  condition  actually  follows  from  the  first  and  (1.5)' 


be  replaced  by  any  larger  number. 

The  solution  of  (l.l)  investigated  by  Atkinson  is  not  the  most  efficient 
multiply  reflected  wave  series  solution  from  the  point  of  view  of  convergence 
properties^  however.   A  more  rapidly  converging  series  solution  of  the  same 
type  existsL J  and  for  this  solution  the  condition  (1.5)  can  be  replaced  by 
the  weaker  condition 

00 

/  |k'(x)/k(x)  jdx  <  00  .  (1.6) 

o 

Moreover^  the  construction  of  this  more  rapidly  converging  series  can 
be  generalized  in  a  straightforward  way  to  the  case  of  a  system  of  first  order 
differential  equations.   An  analogous  series  solution  for  the  system  can  be 
constructed^  and  the  condition  for  the  convergence  of  this  series  is 
correspondingly"  as  weak  as  (1.6)  in  an  appropriate  sense.   This  series  solution 
for  the  system  converges  more  rapidly  and  under  weaker  conditions  than  a  similar 
one  investigated  by  Atkinson L  -I  who  also  generalized  (l.l)  and  (1.2)  to  the 
case  of  a  system  of  ordinary  differential  equations.   A  general  procedure  for 
constructing  a  multiply  reflected  wave  series  solution  for  a  system  of  differential 
equations  is  essentially  contained  in  the  work  of  Keller  and  Keller,!-  -•  of  which 
the  series  in  [6|  is  a  special  case. 

It  is  the  purpose  of  the  present  remarks  to  call  attention  to  the  existence 
of  these  alternative  m\iltiply  reflected  wave  series  solutions  so  that  any  false 
impression  concerning  the  limitations  of  this  type  of  solution  may  be  dispelled. 
Actually,  a  slight  modification  of  the  series  will  provide  a  solution  which  is 

Ft! 

valid  even  when  k(x)  has  finite  jump  discontinuities,  i-  -i 


II.  Some  Definitions  and  Relations 

The  following  are  a  set  of  definitions  and  results  used  by  Atkinson 
and  repeated  here  for  convenience: 


M 


u(x)  =  v(x)  +  w(x)  , 


(2.1) 


with  v(x)  representing  the  total  contribution  of  right-moving  waves  and  v(x) 
representing  the  total  contribution  of  left-moving  waves;  y(x)  is  a  solution 
of  the  differential  equation 


y'  =  (ik  -  |k"^  k')y  , 


(2.2) 


y(x)  =  [k(x)]-i/2 


exp 


X 

i  /  k(x)ds 
<-  o 


(2.3) 


z(x)  is  a  solution  of  the  differential  equation 


lv-1  V.' 


z'  =   (-ik  -  pC     k')z  , 


z(x)   =   [k(x)]   '   exp 


X 

-i  /  k(s)ds 


(2.1+) 
(2.5) 


there  are  also  defined  the  functions 


y(t2.\)  =  y(t2)y"^(\)  , 


-1, 


z(t2,t^)   =  z(t2)z""(t^) 


(2.6) 
(2.7) 


The  fiinctions  v(x)    and  ■w(x)    satisfy  the  integral  equations 


X 


v(x)     =     y(x)   +  I     /     y(x,s)k"^(s)k'(s)tf(s)   ds    ,  (2.8) 


Kx)      =      -  I      /   z(x,s)k"^(s)k'(s)v(s)    ds    .  (2.9) 


They  also  satisfy  the  differential  equations 


v'  =   (ik  -  ■|k"-'-k')v  +  ^k'-Tc'v  ,  (2.10) 


(-ik  -  •|k"\')w  +  |k  \'v  ,  (2.11) 


which  together  with  (2.1)  imply  (l.l). 

III.  The  Multiply  Reflected  Wave  Series  and  Convergence 

The  series  whose  validity  was  investigated  by  Atkinson L  J  is  the 
Neiimann  series  solution  of  (2.8-9).   This  solution  can  also  he  obtained 
by  introducing  a  parameter  X  in  front  of  each  integral  in  (2.8-9);  assuming 
series  expansions  for  v  and  w  in  powers  of  \   and  then  setting  X  equal  to  one. 
The  series  solutions  are 

00 

v(x)  =  y""  u_(x)  ,  (3.1) 


2n 


n=0 


w(x)   =   >    U2^^^(x)  ,  (3.2) 

n=0 


■    -  5  - 
where  the  terms  of  (5.1-2)  satisfy  the  recursion  relations 

OS 

''2n+l^''^   "  -  I  J  z(x,s)k"^(s)k'(s)u2^(s)  ds  ,  (j.j) 

X 

00 

''an^''^  =  I  /y(>^.s)k"^(s)k'(s)u2^_^(s)  ds  (5.4) 


and  u  (x)  =  y(x) .   The  conditions  (1-5-5)  on  the  convergence  of  (5-1-2)  were 


given  by  Atkinson. 


[1] 


In  Section  IV  of  this  paper   a  general  method  for  obtaining  a  multiply 
reflected  wave  series  solution  of  a  system  of  first  order  differential 
equations  will  be  described.   This  method  leads  in  a  natural  way  to  integral 
equations  from  which  the  series  can  be  obtained.   Thus^  if  the  second  order 
differential  equation  (l.l)  is  replaced  by  an  equivalent  pair  of  first  order 
differential  equations 

u'  =   -if  , 
f  =  -ik^u  , 

where  the  fiinctions  u(x)  and  f(x)  are  defined  in  terms  of  v(x)  and  w(x)  by 

u  =  w  +  v  , 

f  =  k(w-  v)  , 

then  the  preliminary  steps  of  the  general  method  described  in  Section  IV  of 
this  paper  will  lead  directly  to  (2,10-11).   The  final  steps  of  that  method 
require  the  definitions 


-  6  - 
v(x)   =  y(x)  p(x)  , 

w(x)   =   z(x)  q(x)  , 

where  v(x)  and  w(x)  are  a  solution  pair  of  (2.10-11).   When  these  are 
substituted  into  (2.10-11)  and  use  is  made  of  (2.2-i^-)  there  results  the 
system  of  differential  equations 

p'  =  ^""Tc'y"  zq.  , 

L  -r  ,  -1  ^ 
q   =  2^  Ti'z  yp  . 

By  integrating  both  sides  of  these  equations  one  can  obtain,  in  particular, 


p(x)   =  c  -  -  /  k"  k'y"  zq  tis  , 


q(x)  =      -  -    J   k"T£'z"  yf   ds  , 

X 


where  c  is  an  arbitrary  constant.   This  leads  at  once  to  a  pair  of  integral 
equations  for  v(x)  and  w(x) : 

00 

y(x)  =  cy(x)  -|  J    y(x,s)k"^(s)k'(s)w(s)  ds  ,    (3-5) 

X 

00 

w(x)  =     -  k         z(x,s)k"^(s)k'(s)v(s)  ds  .         (5-6) 

X 


-  7  - 

Observe  that  w(x)  approaches  zero  as  x  approaches  infinity;  this  same 
condition  is  satisfied  by  w(x)  in  (2.8-9).   If  the  constant  c  is  selected 
appropriately  the  solutions  v(x) ,  w(x)  of  (5.5-6)  can  be  made  identical 
with  the  solutions  v(x) ,  w(x)  of  (2.8-9).  From  a  consideration  of 

_^  Q  [v(x)/y(x)]  =   _^   [v(x)/y(x)]  =  1  it  can  be  seen  that  the  proper 


choice  is 


00 

c  =  1  +  Ij  y"-^(s)k"^(s)k'(s)w(s)  ds  .  (5.7) 


In  order  to  obtain  a  series  solution  of  (5'5-6)  one  can  use  a  Neumann 
expansion  starting  with  v  (x)  =  cy(x) .  The  recursion  formulas  for  the  terms 
of  the  series  are 

CO 

'^2n^''^   =   "  i  J  y(x.s)k"^(s)k'(s)u2^_^(s)  ds  ,         (5.8) 


00 

''2n+l^''^  ^      '  k    J   2(x.s)k"^(s)k'(s)u2^(s)  ds  .  (3-9) 


X 


The  series  solutions  of  (5-5-6)  are  then 


00 


v(x)  =   >   ^sn^""^  '  (^'^^^ 

n=0 

00 

w(x)  =  /,  ^2n+l^''^  •  (5.11) 

n=0 

It  can  be  shown  by  means  of  a  standard  procedure  L  -1  that  the  Neiomann 
series  for  (3«5-6)  converges  rapidly  under  the  condition  (I.6),  for  the  integral 


-  8  - 

equations  (3.5-6)  are  of  the  Volterra  type.   It  is  possible  also  to  prove  this 
convergence  property  by  means  of  a  majorizing  series  of  the  type  used  by 
Atkinson. L  J 

Following  Atkinson  one  can  define 

P(x)   =  y"  (x)v(x),     q(x)   =   z"  (x)w(x) 
so  that  the  integral  equations  (5- 5-6)  become 

00 

p(x)  =  c  -  I  J  y"V\'zq  ds  (5.12) 


X 


q(x)   =   -  I  /  z"V\'yp  ds  ,  (3.13) 


X 


and  the  recursion  relations  (3.8-9)  become 

00 

P2,W  =  -i  /y-V\'zp2^_^ds  ,  (5.11.) 


X 


P2,,i(x)   =  -i  /z-V^k-yp^^cis  .  (3.15) 


X 


The  solutions  p  (x)  of  (5.1^1-15)  are  terms  in  the  series 


p(x)   =  2_j''2n^'')  '  ^^-^^^ 

n=0 


.(x)   =  ^P2n.lW  •  (5.17) 

n=0 


-  9  - 

The  convergence  of  the  series  (5.I6-I7)  depends  upon  the  convergence  of 

00 

n=0 

in  the  unit  circle  in  the  complex  X   plane.  As  in  Atkinson's  analysis  this 
series  is  majorized  by 


r(x)  =  Y^  ajx)x''   ,  (3.18) 


n=0 


where  a   (x)  =  c  and  the  remaining  a   (x)  are  defined  by 


snW  - 1  r  w'^^'\  s._-,<is  >       (5.19) 


2n'  '     2    J       '  '   2n-l 

X 


cr^^^lCx)  =  I  /   |k-\'|  a^^ds  .  (5.20) 


The  absolute  convergence  of  (3.I8)  is  equivalent  to  that  of 


n=0 


The  series  (3.21)  can  be  summed  through  a  consideration  of  the  differential 
ecjiation  which  it  satisfies  and  the  conditions 


0(ooj  =      c    , 


^^■^  ^  |k-\' 


=   0 


-  10  - 
The  sura  is  found  to  be 

o-(x)   =  c  cosh  I  X   /   |k"\'|ds|  .  (3-22) 

Since  this  function  is  an  entire  function  of  X  the  majorizing  series 
converges  for  all  X  and  in  particular  for  X=  1.   Thus^  the  series  (5-16-17) 
converges  under  the  condition  (1.6). 

* 
IV.  The  Case  of  a  System  of  Differential  Equations 

The  differential  equations  (2.10-11)  for  a  propagating  wave  can  be 

generalized  to  a  system  of  differential  equations  which  can  then  be  replaced 

by  a  system  of  Volterra  integral  equations.   Let  A(x)  be  an  n  by  n  matrix 

and  u(x)  an  n  dimensional  column  vector.   Consider  the  differential  equation 

u'(x)  =  A(x)u(x)  .  (^.1) 

Let  K(x)  be  an  n  by  n  matrix  which  is  nonsingular  for  almost  all  values 
of  X  such  that  K  AK  is  diagonal  with  diagonal  elements  k  (x) .   After  the 
substitution 

u(x)   =  K(x)v(x) 
there  results  in  place  of  (i|-.l)  the  differential  equation 


v'(x)  =  [-K"-'-(x)K'(x)  +  k"^(x)A(x)k(x)]v(x)   .    {h.2) 


The  solution  given  in  this  section  is  due  to  Keller  and  Keller. 


p] 


11  - 


The  differential  equation  (i4-.2)  is  a  generalization  of  the  system  (2.10-11). 

Let  the  diagonal  elements  of  the  matrix  k'  (x)K'(x)  be  7  (x) .  Then 
define  the  matrix  Y(x)  as  the  solution  of  the  differential  equation 

Y'(x)   =   [-r(x)  +  K-^(x)A(x)K(x)]Y(x)  ,  (U.3) 

where  r(x)  is  the  diagonal  part  of  k"  (x)K'(x).  The  matrix  Y(x)  is  diagonal 
and  can  be  calculated  explicitly  since  the  matrix  factor  of  Y(x)  on  the  right 
of  (4.3)  is  diagonal.  The  solution  Y(x)  is  determined  except  for  an  arbitrary 
constant  diagonal  matrix  factor  which  can  be  taken  to  be  the  identity.  The 
n-th  diagonal  element  of  Y(x)  will  have  the  form 


Y  (x)  =  exp 
n 


X 

/  7^(s)ds  +  /  K^(s)  ds 


■-   o 


(h.k) 


Now  define  the  vector  p(x)  by 


v(x)   =  Y(x)p(x) 

and  substitute  into  (i<-.2).  There  resiolts,  with  the  use  of  (^.3)^  the  differential 
equation 


p'(x)   =  Y"^(x)H(x)Y(x)p(x) 


(i+.5) 


where  H(x)  is  K  (x)K'(x)  with  each  of  its  diagonal  elements  replaced  by  zero. 
An  integral  equation  for  p(x)  can  be  obtained  from  (k.^)   by  integrating  both 
sides 


00 

p(x)  =  c  -  /  Y"-^(s)H(s)Y(s)ds  , 


{h.6) 


-  12  - 

where  c  is  an  arbitrary  constant  vector. 

Left-moving  and  right-moving  wave  components  of  the  wave  vector  v(x) 
can  be  defined  according  to  whether  the  imaginary  part  of  the  corresponding 

X 

quantity  /  k    (x)ds  in  (4A)  is  a  decreasing  or  an  increasing  function  of  x. 

o 
To  ensure  that  no  left -moving  wave  components  are  present  at  x  =  ^o  one  can 

choose  the  constant  vector  c  in  (4.6)  so  that  its  components  corresponding 

to  the  left-moving  wave  components  of  v(x)  are  zero. 

The  integral  equation  (h.S)    is  of  the  Volterra  type;  therefore^  its 
Neumann  series  expansion  converges  rapidly  and  under  relatively  weak 
conditions.   If  the  elements  of  a  matrix  M  are  M. .,  a  convenient  norm  for 
the  matrix  is 

-  I1/2 


M 


Im.  . 


(4.7) 


For  the  series  solution  of  (4.6)  the  condition  for  convergence  which 
corresponds  to  (I.6)  is  [6]: 


Y""'"(s)h(s)y(s)  II  ds  <  ~.  (4.8) 


When  the  differential  equation  (4.1)  describes  an  energy  conserving  propagating 

wave  the  diagonal  matrix  Y(s)  is  unitary,  and  in  this  case  the  condition  (4.8) 

becomes 

00 

f     II  H(x)||  ds  <  ~  .  (4.9) 


-  15  - 

It  is  interesting  to  note  that  if  the  frequency  w  is  displayed  by 
replacing  the  matrix  A(x)  in  (i+.l)  with  (jjA(x)  ,  the  error  estimate  of  any- 
partial  sum  of  the  series  solution  of  (k.l)   based  on  the  Neumann  series 
solution  of  (i|-.6)  depends  only  on  the  quantity 

00 

/||H(s)||ds 
o 

and  thus  is  independent  of  cj.  The  implication  of  this  fact  is  that  the 

series  solution  will  converge  uniformly  with  respect  to  u. 
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